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Many organisms live in confined fluidic environments such as the thin liquid layers on the skin of host organisms or in partially- 
saturated soil. We investigate the collective behaviour of nematodes in a thin liquid layer, which was first observed by Gray and 
Lissmann, [J. Exp. Biol. 41, 135 (1964)]. We show experimentally that nematodes confined by a thin liquid film come into 
contact and only separate again after some intervention. We attribute this collective motion to an attractive force between them 
arising from the surface tension of the layer and show that for nearby nematodes this force is typically stronger than the force 
that may be exerted by the nematodes' muscles. We believe this to be the first demonstration of the "Cheerios effect" acting 
on a living organism. However, we find that being grouped together does not significantly alter the body stroke and kinematic 
performance of the nematode: there are no statistically significant changes of the Strouhal number and the ratio of amplitude to 
wavelength when aggregated. This result implies that nematodes gain neither a mechanical advantage nor disadvantage by being 
grouped together; the capillary force merely draws and keeps them together. 
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1 Introduction 

Organisms need to move for a variety of reasons including to 
avoid predators, to find nutrients or to find a mate. This mo- 
tion often takes place in fluid environments and hence may in- 
volve hydrodynamic interactions, which have been studied ex- 
tensively in recent years. At the microscopic scale, hydrody- 
namic interactions between or gani sms can result in the genera- 
tion of large-scale circulationst^Zl, a highly concentrated layer 
in a shear field ^ or the non-Brownian motion of a test parti- 
cleK. This collective behaviour can also enhance mixing and 
transport of chemical nutrientsEEl. In the majority of recent 
studies, collective movements take place in the bulk of flu- 
ids and hence owe their origin to bulk hydrodynamic interac- 
tions. Here we report the temporary aggregation of nematodes 
caused by the presence of a liquid-gas interface. Ostensibly, 
this is similar to the dancing of Volvox (algal colonies) at the 
surface of water ^, but the physics is very different. 

Nematodes locomote by creating undulatory movement 
similar to other invertebrates. The propagating body- wave is 
generated by contracting and relaxing muscles under the con- 
trol of the neuromuscular system. The maximum force that 
nematodes are capable of generating has been measured in 
micro-fabricated devices, and shown to be on the order of 1 
This limiting force suggests that the application of a 
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force with magnitude greater than 1 /iN may change or con- 
strain the organism's behaviour. 

It is often observed that nematodes migrate up the side 
of their culture jars and move into a thin liquid film region. 
When this occurs they are also often observed to aggregate 
together Qualitatively, such collective behaviour has been 
found in many types of nematodes (such as the hookworm 
larvae and parasitic nematodes ^ ^ ) and has been shown not to 
originate in chemical interactions^^. In this article we shall 
suggest that this aggregation is caused by lateral capillary 
forces — similar to the "the Cheerios effect" observed with 
floating objects. 

Capillary forces arise from the interfacial tension of fluid 
phase boundaries, and act along the interfaces. They can, 
therefore, be important for small organisms living close to 
fluid interfaces; they allow, for example, the pond skater to 
walk on the surface of water When a liquid film is thinner 
than the diameter of a nematode body the nematode experi- 
ences both vertical and lateral forces because of the surface 
tension of the interface. The vertical component of surface 
tension pulls its body down against the substrate, lowering the 
lateral slipping during locomotion and allowing them to crawl 
on a surface. The horizontal component of surface tension, 
meanwhile, will vanish when the nematode is alone but may, 
in principle, lead to an interaction with a nearby nematode. 

In this paper, we shall investigate the collective behaviour of 
nematodes in a thin liquid layer. We characterize the locomo- 
tory performance of nematodes by measuring the geometric 
parameters of their body stroke and forward moving veloc- 
ity. We compare this locomotory performance for the cases of 
aggregated and individual nematodes. We then present a sim- 
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Fig. 1 Two nematodes (Panagrellus redivivus) merging: (a) two long nematodes and (b) two small nematodes. Typically, it takes about 1 sec 
for them to merge completely from the initial point of touching. 




Fig. 2 Various stages of aggregating, joint, and separating 
nematodes. Two nematodes are aggregated at ^ = 6.7 sec, and stay 
together for more than 10 sec. However, they are hit by a third 
nematode at ? = 17.7 sec ultimately leading to the separation of all 
three nematodes. (Supplementary Video 1$) 

plified theoretical model to estimate the typical lateral forces 
due to surface tension and compare it with the typical forces 
generated by nematodes during locomotion. Our study ratio- 
nalizes the physical basis of the aggregation but demonstrates 
ultimately that it changes little about the locomotion stroke. 

2 Experimental methods 

2.1 Nematodes crawling on a thin fluid 

Panagrellus redivivus was used as the model organism in our 
experiments. P. redivivus, often called the micro- worm, is a 
free living, bacteriophagous nematode. It can be found in soil, 
vegetation, and dung^^, demonstrating its ability to move in 
various fluidic conditions. Cultures of P. redivivus were main- 
tained at a temperature of 18 °C in a potato medium. P. redi- 
vivus grows to around 1 mm in length L and 50 /im in diameter 
2rc. 

Nematode cultures were placed in sealed glass chambers. 



Over several hours some nematodes migrate up the side of the 
glass. We characterize the motion of nematodes a few cen- 
timetres above the bulk of the culture medium. At this height, 
only a thin liquid layer remains on the side glass wall, which 
is quite different from the slurry bulk culture below. Videos 
of nematodes are taken using a Sony Handycam with a macro 
lens and Virtualdub software was used to format video. Data 
are taken from the same chamber with the same ambient air 
temperature to avoid any thermotatic behaviours. The MAT- 
LAB image processing toolbox is used to analyze footage and 
process data. Peak-to-peak amplitude A and wavelength A of 
the nematodes are measured at each stage of the aggregation 
process. Nematode speed V is recorded by tracking the cen- 
troid of the body typically over a few periods of undulation. 
The frequency / of nematode undulation and the length L of 
each nematode are also measured. 

Figure [T] shows two nematodes coming together as they 
crawl in proximity to each other. Typically, their heads touch 
and their forward motions draw the rest of the bodies together. 
Once they are grouped, the capillary force makes it difficult 
for the nematodes to separate themselves. We observe that ne- 
matodes stay together from as little as a few seconds to as long 
as a minute depending on whether they are hit by others (see 

We often observed events in which two nematodes of dif- 
ferent sizes aggregate as in Fig. [TJb). Before two nema- 
todes merge, the individuals have different body strokes and 
speed. Typically, smaller nematodes have shorter wavelength 
and amplitude and lower speed. However, after they aggre- 
gate they form a single body form — the smaller nematode 
adopting the body form of the larger nematode. 

To characterize the aggregation behaviour, the peak-to-peak 
amplitude over wavelength and non-dimensional Strouhal 
number (the ratio of lateral to forward velocities) were calcu- 
lated for individual and paired nematodes. The ratio of wave- 
form amplitude to wavelength is a primary indicator of chang- 
ing body stroke in swimming organisms while the Strouhal 
number gives an indication of the efficiency of the locomo- 
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Fig. 3 (a) Amplitude versus wavelength, and (b) /A versus forward 
velocity, V, for individual and grouped nematodes. In each case, the 
lines correspond to the mean ratios for individual (dashed line) and 
grouped (solid line) nematodes. 



tion strokJ^. Data are taken from 24 individual nematodes, 
and 23 cases of grouped nematodes and the 95% confidence 
intervals for mean ratio and mean Strouhal number are con- 
structed for each condition (individual or grouped). The 95% 
confidence interval for the mean ratio of amplitude to wave- 
length was (0.79,0.95) amongst individuals whilst amongst 
grouped nematodes the corresponding confidence interval was 
(0.76, 0.90). Since these confidence intervals overlap, we con- 
clude that there is no significant difference between the mean 
ratio of stroke amplitude to wavelength. However, we note 
that these values are higher than the value reported in crawl- 
ing nematodes on a soft agar gel (0.63±0.03)IS1. The 95% 
confidence interval of mean Strouhal number for individu- 
als was (0.54,0.76) while that for grouped nematodes was 
(0.57,0.84). Since these confidence intervals overlap, we 
again conclude that there is no statistically significant differ- 
ence between the mean Strouhal number for grouped and in- 
dividual nematodes. In conclusion, the features of their loco- 
motion are not changed by being aggregated. Figure |3] demon- 
strates this similarity in the locmotion characteristics for indi- 
vidual and grouped nematodes graphically. 

To characterize the dominant mechanism for nematode ag- 
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Fig. 4 Collision time vs. crawling time. Collision time is measured 
as the time interval from point- wise touch to grouped two-body and 
crawling time is defined as the ratio of body length to crawling 
velocity. 



gregation, we consider two possible time-scales associated 
with the aggregation. One is the capillary time-scale account- 
ing for the surface tension force that draws them together, de- 
fined as Tcapiiiary = 6;r/irc/7 where /i is the fluid viscosity (4.0 
cP for a filtered potato medium) and 7 is the surface tension 
coefficient (note that this holds even for zipping bodies ^^). 
The other is the crawling time defined as Tcrawl = L/V. 

Experimentally, we define the collision time as the time in- 
terval from the initial contact to the grouped two-body. The 
grouped two-body is referred as a configuration when the two 
nematodes are touching along their body length but do not 
necessarily have their heads and tails in the same phase. For 
example, the second column in Fig. [TJa) is when two nema- 
todes make the initial contact and the sixth column represents 
the grouped two-body. Hence, in Fig. [TJa) the collision time is 
1.88 sec. A comparison between crawling and observed col- 
lision time scales is shown in Fig. |4] We observe that the 
capillary time scale is expected to be on the order of 10~^ 
sec and hence is significantly smaller than the measured col- 
lision time scales. Furthermore, the scaling for Tcapiiiary does 
not explain the observed relationship between size, speed and 
collision time. However, the crawling time depends on the 
body size and the crawling velocity, and is on the order of the 
measured collision times. From this, it seems clear that the 
initial aggregation is driven by crawling, rather than by cap- 
illary forces. We therefore hypothesize that aggregation ini- 
tially occurs because of collisions but that, once aggregated, 
the capillary force keeps the two nematodes in one body form. 
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3 Model 

We now seek to rationalize the observation that nematodes 
seem to be bound together after having come close to one an- 
other. We begin by considering the relative importance of the 
various physical forces in the system. Firstly, we note that 
the Bond number, Bo = pgl} jy 0(0. 1): thus surface ten- 
sion dominates gravity for nematodes. In what follows, there- 
fore, we shall neglect gravity entirely. We also note that the 
Reynolds number, which measures the relative importance of 
inertia to viscous effects. Re = pVL/ji ~ (9(10~^) and so liq- 
uid inertia is entirely negligible here. The insignificance of 
inertia means that velocities are proportional to forces and it is 
enough to calculate the typical interaction force between two 
nearby nematodes. In principle, it would then be possible to 
calculate the perturbation caused to their motion by this force, 
along the lines of capillary zippering^^, though we have al- 
ready seen that this perturbation must be small. 

To estimate the typical interaction force between two nema- 
todes we calculate the typical force that 'immersion forces' 
can exert on two horizontal cylinders. This has not been 
studied before although the force between floating cylinders 
has been calculate d as has the force between immersed 
spheres (The case of immersed objects is different from 
that of floating objects because in the former case the substrate 
provides whatever vertical force is required to satisfy the ver- 
tical force balance condition. The vertical force balance must 
be considered separately for floating objects.) We consider the 
idealised setup illustrated in Fig. [5] (a), which is qualitatively 
similar to that observed for an actively moving nematode (see 
Fig.|5](b)). 




Fig. 5 (a) Schematic of cross-sectional view of two worms, (b) 
Experimental observations of meniscus around a moving nematode. 
(Supplementary Video 2$) 



Both between the cylinders and outside, the angle of the 
meniscus to the horizontal satisfies (j) = 71 — {6c -\- y/) at the 
contact line, where 6c is the contact angle on the cylinder. El- 
ementary geometry then shows that the radii of curvature of 
the menisci in the inside and outsideregions are 
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The total horizontal force consists of two components: a force 
due to surface tension acting at the contact line and a pressure 
force acting along the wetted perimeter. The horizontal force 
per unit length due to surface tension acting at the contact line 
is given by 



Fy = 7(cos 0; — cos 



(2) 



Here, a positive force corresponds to an attractive interaction. 
There is also a contribution to the force from the difference in 
pressures between the inner and outer menisci. The pressure 
in the fluid in these regions is given by p — Pa 
inward pressure force per unit length is then 
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Therefore, the total horizontal force per unit length on the 
cylinder is 

^ =cos0—cos0o + ^(l-cosVA/)-^(l-cosv/-o). (4) 
7 Ri Ro 

In eqn. ^ we have the force for given values of W and 
However, we need to estimate these angles to estimate the total 
horizontal force. To obtain an equation relating i//^ and Wi we 
shall assume equal capillary pressures in the inside and outside 
regions, since then Ri = Ro. This may be justified physically 
since it is unlikely that nematodes will be in contact with the 
solid substrate all along their length. Hence any difference in 
capillary pressure would be quickly equalized. We shall also 
make the simplifying assumption that the nematodes and the 
glass that they sit on are both perfectly wetting so that 6c = 0. 

With these simplifications, we find that the nonlinear terms 
cancel and we have simply 



cos Wo 



: sin Xf/i — 1 



(5) 



where 5 = d/2rc > 1. Therefore, the total horizontal force per 
unit length simplifies to 



Ftot 
7 



: sin Xj/i — cos Wi — 1 



Ri 



(6) 
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Fig. 6 Force-displacement curves for different values of i///: i/// = 1 
(red), \j/i = \ .5 (green) and — 2 (blue). Here 3 = d/2rc and we 
have assumed that = 0, as is likely to be the case for the worms. 



nematodes. This result indicates that the body form is not sig- 
nificantly changed by the capillary force and, further, that no 
mechanical advantage is achieved. Nevertheless, grouped ne- 
matodes are observed to stay together for a long time (more 
than a few minutes in some cases). It is known that many 
other organisms move and live in circumstances where capil- 
lary forces are important and so it seems likely that, in some 
of these instances at least, the capillary force may be exploited 
to the organism's benefit, as has already been shown for indi- 
vidual water treaders, Mesovelid^. 

5 Acknowledgements 

D.V. is supported by an Oppenheimer Early Career Fellow- 
ship. 



For this force to be attractive, we require that | sin i/// — 
cos i/A; > 1, that is, sin(i/A; - a) > where tana = 5/2. 

For a given displacement, i.e. given 5, we have that the force 
is attractive provided that 

f 5 \ 5 

For simplicity, we shall assume that the contact line is 
pinned so that i//^ is constant as the cylinders move closer to- 
gether. By plotting the force-displacement curves we find that, 
as expected, the attractive force increases as the two worms 
get closer together. Some typical results are shown in Fig. 
[6] These plots also demonstrate that the force only acts over a 
lengthscale comparable to the radius of the worms; the force is 
relatively short ranged. This is consistent with the observation 
that zipping is driven by collision, rather than being driven by 
surface tension. Finally, the total horizontal force on a worm 
of length 1 mm could easily reach the order of 50 /iN, and 
consequently overwhelm the typical muscle force that can be 
exerted by a single nematode ~ 0(1) /iN. 

4 Discussion 

We have studied the collective behaviour of nematodes in thin 
liquid layers, thereby uncovering the possibility of interactions 
between locomoting organisms without hydrodynamic effects. 
This effect is observed even when the nematodes are of very 
different sizes. We have shown that this interaction is not ini- 
tially driven by surface tension forces but that once an inter- 
action happens surface tension forces ensure that the nema- 
todes remain aggregated. However, we have also shown that 
the characteristics of locomotion (namely the Strouhal number 
and body form) do not differ between individual and grouped 
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